DECISION FEEDBACK EQUALIZERS WITH CONSTRAINED FEEDBACK 
TAPS FOR REDUCED ERROR PROPAGATION 



Related Applications 

This application claims the benefit of U.S. 
Provisional Application 60/457,141 filed in March 24, 
2003. 

Technical Field of the Invention 

The present invention relates to decision 
feedback equalizers having constrained taps. 

Background of the Invention 

The presence of multipath interference in a 
broadcast channel, such as a broadcast television 
channel, has been a persistent obstacle to good reception 
by a receiver. The broadcast television channel is a 
relatively severe multipath environment due to a variety 
of conditions that are encountered in the channel and at 
the receiver. Channels are characterized by channel 
impulse responses which may be several hundreds of 
symbols long, so that strong interfering signals may 
arrive at the receiver both before and after the largest 
amplitude signal. In addition, the signal transmitted 
through the channel is subject to time varying channel 
conditions due to the movement of the transmitter and/or 
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due to signal reflectors, airplane flutter, and, for 
indoor reception, people walking around the room. If 
mobile reception is desired, movement of the receiver 
must also be considered. As is known, intersymbol 
5 interference (ISI) is also an obstacle to good reception. 
Multipath and intersymbol interference 
adversely affects the ability of the receiver to 
correctly receive the symbols transmitted by the 
transmitter. Therefore, equalizers are used in receivers 
10 in order to cancel the effects of multipath and 

intersymbol interference and thereby improve signal 
reception . 

A decision feedback equalizer (DFE) is an 
attractive type of equalizer because it theoretically 

15 performs better than a linear equalizer. The decision 
feedback equalizer typically includes a feed forward 
filter that includes a plurality of feed forward taps, a 
feedback filter that includes a plurality of feedback 
taps, some type of decision device such as a slicer or 

20 trellis decoder, a summer, and a tap weight calculator 

that sets the tap weights for the feed forward filter and 
the feedback filter. 
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The received signal is provided to the feed 
forward filter. The outputs of the feed forward filter 
and the feedback filter are summed by the summer. The 
output of the summer forms the output of the decision 
5 feedback equalizer. The decision device decides correct 
values for the output symbols of the summer and feeds 
these decisions to the feedback filter. The tap weight 
calculator estimates the channel impulse response based 
on the received signal and the decisions from the 

10 decision device. The tap weight calculator then 

determines the optimal (MMSE - minimum mean squared 
error) tap weights based on this channel impulse response 
estimate and the estimated signal-to-noise ratio that 
characterizes the received signal, and adjusts the feed 

15 forward taps of the feed forward filter and the feedback 
taps of the feedback filter according to these tap 
weights . 

If the tap weights of the feedback filter of 
the decision feedback equalizer are relatively large, and 
20 if the decision device of the decision feedback equalizer 
makes an error, a problem known as error propagation may 
arise. Error propagation can significantly reduce the 
performance of a decision feedback equalizer. In cases 
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of severe intersymbol interference (ISI), this reduction 
in performance can be so large that steps must be taken 
to mitigate the error propagation. 

In order to reduce error propagation, it is 
known to calculate the MMSE tap weights for the decision 
feedback equalizer subject to a constraint on the tap 
weights of the feedback filter, thereby limiting the tap 
weights of the taps of the feedback filter to values that 
reduce error propagation. In order to calculate such 
constraint limited tap weights, it is typically necessary 
to differentiate the constraining function. However, 
certain desirable constraint functions, such as the 1- 
norm (that maintains the sum of the absolute values of 
the feedback tap weights below a selected value) are not 
dif f erentiable and, therefore, cannot be used for this 
purpose . 

The present invention overcomes one or more of 
these or other problems . 

Summary of the Invention 

According to one aspect of the present 
invention, a method of supplying tap weights to taps of a 
decision feedback equalizer comprises the following: 
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determining a channel impulse response of a channel 
through which the decision feedback equalizer receives a 
signal; determining constrained tap weights based on the 
channel impulse response and a dif f erentiable tap weight 
5 constraint function, wherein the dif f erentiable tap 

weight constraint function is an approximation of a non- 
diff erentiable tap weight constraint function; and, 
supplying the constrained tap weights to the decision 
feedback equalizer . 

10 According to another aspect of the present 

invention, a decision feedback equalizer comprises a feed 
forward filter, a decision device, a feedback filter, a 
summer, and a tap weight determiner. The feed forward 
filter applies feed forward taps to a signal to be 

15 equalized. The feedback filter applies feedback taps to 
an output of the decision device. The summer sums 
outputs from the feed forward filter and the feedback 
filter to provide the output of the decision feedback 
equalizer, and the output of the decision feedback 

20 equalizer is provided to the decision device. The tap 

weight determiner determines constrained tap weights and 
unconstrained tap weights based on a minimum mean squared 
error between the signal to be equalized and the output 
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of the decision device. The constrained tap weights are 
determined by constraining the minimum mean squared error 
according to a tap weight constraint function, and the 
tap weight constraint function is dif f erentiable and is 
5 an approximation of a non-diff erentiable tap weight 
constraint function. The tap determiner supplies the 
constrained tap weights to the taps of the feedback 
filter and the unconstrained tap weights to the taps of 
the feed forward filter. 

10 According to still another aspect of the 

present invention, a method of supplying tap weights to 
taps of a decision feedback equalizer comprises the 
following: determining a channel impulse response of a 
channel through which the decision feedback equalizer 

15 receives a signal; determining a constraint value M as a 
function of an output of the decision feedback equalizer; 
determining constrained tap weights based on the channel 
impulse response and a tap weight constraint function 
having the constraint value M; and, supplying the 

20 constrained tap weights to the decision feedback 
equalizer . 

Brief Description of the Drawings 
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These and other features and advantages will 
become more apparent from a detailed consideration of the 
invention when taken in conjunction with the drawing in 
which : 

5 Figure 1 shows a baseband equivalent discrete- 

time equalizer system useful in describing the present 
invention; 

Figure 2 is a graph illustrating differences 
between several modified constraints on tap weight values 
10 for the feedback taps of a decision feedback equalizer 
and the unmodified constraint; 

Figure 3 illustrates a decision feedback 
equalizer according to the present invention; 

Figure 4 shows an apparatus for adaptively 
15 determining the constraint value M; 

Figure 5 is a graph that illustrates the 
adaptation of the constraint value M as a function of the 
mean squared error at the output of the decision feedback 
equalizer; and, 
20 Figure 6 is a graph that illustrates an 

alternative adaptation of the constraint value M as a 
function of the mean squared error at the output of the 
decision feedback equalizer. 
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Detailed Description 

A baseband equivalent discrete-time equalizer 
system 10 is shown in Figure 1. The sequence I[k] 
5 represents the transmitted symbols which are assumed to 
be independent and equally probable with zero mean and a 

variance £ T = E|l[£:]| 2 j, where E is the expectation 

operator. For 8-VSB such as used in digital television 
according to the 1996 ATSC digital television (DTV) 

10 standard, the transmitted symbols are taken from the set 
A = {±1, ±3, ±5, ±7} . 

In continuous time, the symbols I[k] are 
transmitted through the multipath channel 12. The 
multipath channel 12 is denoted c[k] and contains M paths 

15 with delays r ± and complex gains a ± . Including the 

effect of pulse shaping at the transmitter, the channel 
impulse response of the channel 12 is given by the 
following equation : 

M 

20 c(t) = £ a ± q(t - r ± ) . (1) 
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The path considered to be the "main" path is always 
assigned a delay r of zero. In an ATSC 8-VSB system, the 
transmitted pulse q(t) has a square-root raised cosine 
spectrum with a nominal bandwidth of one-half the symbol 
5 rate, but is shifted up in frequency by one-fourth the 

symbol rate. Therefore, the transmitted pulse shape q(t) 
is complex. 

The discrete-time channel c[k] is found by 
sampling c(t) at the symbol rate. A pulse shaping filter 

10 14 at the receiver is denoted q[k) and is obtained by 

sampling q(t) at the symbol rate. Because q* (— t) = q(t), 
the pulse shaping filter 14 is matched to the transmitter 
filter (where * denotes complex conjugation) . It is 
assumed that q[k] has finite length and is given in 

15 vector form by q = [q[— L q ] ... q[0] ... q[L q ]] T . The 

composite channel is a combination of the channel and the 
pulse shaping and is defined as h[k} = c[k] * q(k] (where 
* indicates convolution) . It is assumed that the 
composite channel also has a finite length and is given 

20 in vector form by h = [h[— L ha ] ... h[L hc ]] T . The length of 
the composite channel is L h + 1, where L h = L ha + Lh C / 
where L ha is the number of taps prior to the main path, 
and where L hc is the number of taps after the main path. 
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The noise sequence rj[k] is introduced as shown in Figure 1 
and is complex additive white Gaussian noise (AWGN) with 
a power spectral density of No- 

The received signal z[k] is given by the 
5 following equation: 

z[k] = h[k] * I[k] + v [k], (2) 

where v [k] = rf [k] * q[k] is colored complex noise. For ease 

10 of explanation, it is assumed that the equalizer is a 
real equalizer so that the real part of the signal is 
taken prior to a feedforward filter 18 of a decision 
feedback equalizer 20. A subscript R is used herein to 
indicate the real part of a given quantity. Therefore, 

15 the input signal to the decision feedback equalizer 20 
may be written according to the following equation: 

z R [k] = Re{z[k] } = h R [k] * I[k] + v R [k], (3) 

20 where h R [k] = Re{h[k)} and v R [k] = Re {v [k] } . Because the 

feed forward filter 18 of the decision feedback equalizer 
20 contains (N FF + 1) taps, (N F f + D consecutive samples 
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of z R [k) are assembled into a vector according to the 
following equation: 



z R [k] = H R l[k] + v R [k] 



(4) 



10 



where z R [k] = [z R [k] ... z R [k — N F f]] t is the vector at the 
input of the decision feedback equalizer 20, I[k] = 
[I[k + L ha ] •-• I[k - L hc _ N FF ]] T is a vector of 
transmitted symbols, and V R [ic] = [v R [k] ■•• v R [k - N FF ] ] T is the 
real part of the colored complex noise. It is known that 
v R [k] is white. The channel matrix may be given by the 
following equation: 



15 



0 h 



0 
0 



(5) 



The colored complex noise vector is given by 
the following equation: 



20 



v[k] = Qr\[k] 



(6) 
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where n W = [?J [k + L q ] ••• 77 [fc - L g 



FF 



f and Q 



is given by 



the following equation: 



Q = 




T 

q 



(7) 



T 

q 



5 



The remainder of the decision feedback 



equalizer 20 includes a decision device 22, a feedback 
filter 24, and a summer 26. The summer 26 adds the 
output of the feed forward filter 18 and the output of 

10 the feedback filter 24, and the decision device 22 

processes the output of the summer 26 by determining 
correct values for the symbols at the output of the 
summer 26. The feed forward filter 18 applies feed 
forward taps to z[k]. The feedback filter 24 applies 

15 feedback taps to the symbol decisions from the decision 
device 24. The output of the summer 26 is the output of 
the decision feedback equalizer 20. 



and the feedback filter 24 has N FB taps. The coefficient 
20 vectors for the feed forward filter 18 and the feedback 



The feed forward filter 18 has (N FF + 1) taps 
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filter 24 are g F = [g F [0] g F [N FF ]] T and 

9b = [^st 1 ! Qb^ n fbW * respectively. At the output of the 
decision feedback equalizer 20, the symbol estimate is 
given by the following equation: 

5 

I[k - K] = g T F z R [k] + g T B I[k - K - 1], (8) 



where l[k] = ••• l[k + 1 - N FB ]\ T are the decisions from N FB 

previous symbols. It is noted that T[k] <= A. The integer 
10 K is defined to be the cursor, which identifies the 

symbol being estimated in I[k]. The decision device 22 
may be either a standard 8-level slicer (quantizer) or a 
trellis decoder as described herein. 

The following equations may be assumed: 



15 



h I 

R ZR x B = H*A K (11) 
r IZR = K R t>K + L ha d2) 
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where is the inverse of the signal-to-noise ratio. 

* i 

The vector 5^ contains all zeros except for a one in 
component n + 1. When N FB < L hc + N FF - K, then 



A* = 



0 



{L ha + K + l)xN FB 
I N FB 



o 



(N FF + L hc -K- N FB )xN FB 



(13) 



10 



In this system of equations, is a zero matrix with m 

rows and n columns, and I n is the identity matrix. 

The mean square error (MSE) at the equalizer 
output is given by the following equation: 



MSE(g F , g B ) = £pl[* — K] — l[k - K]f^ 



= Z x ~ ^xz R + 2g£R 2iRjB g B + 9 t f R ZrZr 9f + 4i9 t b 9b (14) 



15 where E is the expectation operator. The minimum mean 
square error decision feedback equalizer tap weights 
(g"F,gB) are calculated by differentiating equation (14) 
with respect to g F and g B , setting the respective 
derivatives equal to zero, and then simultaneously 
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solving the two equations for g F and g B . This method is 
well known. 

In considering error propagation, the error at 
the output of the decision feedback equalizer 20 is given 
5 by e[k] = I[k] - I[k] . Substituting equations (4) and (8) 
into e[k] yields the following equation: 

e[k] = g£H^I[Jc + K] + g T F v R [k + K] + g T B l[k - 1] - I[Jc] (15) 

10 The channel matrix can be partitioned in accordance with 
the following equation: 

H i? = t H *F h R t] H RP1 U RP2 ] (16) 

15 where H RF has K columns, h ( / t) is a single column (and is a 
time-reversed and possibly truncated channel vector) , H RP i 
has N FB columns, and H RP2 has L h + N FF - K - N FB columns. 
Also, I[k + K] can be similarly partitioned according to 
the following equation: 

20 

I[k + K] = [l F [kf l[k] l P1 [kf I P2 [kf] (17) 
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where similarly I F has K columns, I T is a single column, 

\ T P1 has N FB columns, and 1 T P2 has L h + N F f - K - N FB columns. 
According to equations (15) -(17), the equalizer error may 
be given by the following equation: 

5 

e[k] = 9 T F H RF I F [k] + (g^h ( / t} - l)[[k] + 9 T F H RP1 I P1 [k] 

+ 9%H R p21p2M + 9 T F v R [k] + g T B l[k - 1] . (18) 

Because R Zr i b = h kpi and Y iz R = h if ^ ' the following 
10 equation is obtained: 

g B = "H RP1 r g F (19) 

(the feedback filter taps g B results from a convolution of 
15 the channel and the feed forward filter taps g F ) and the 
equalizer error is now given by the following equation: 

e[k] = g T B (l[k - 1] - I P1 [*]) + (g T F h% t] - l)l[k] 

+ g T F K RF l F [k] + g T F n RP2 I P2 [k] + g T F v R [k] (20) 

20 
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The first term of equation (20) is due to decision 
errors, the second term of equation (20) is due to 
equalizer bias, the third and fourth terms of equation 
(20) are residual intersymbol interference (ISI), and the 
5 last term of equation (20) originates from the noise. 

Equation (20) demonstrates how error 
propagation occurs in the unconstrained equalizer. If a 
single nearest-neighbor decision error (magnitude 2) 

exists in I[k - 1] , and when that error reaches the 
10 position of the maximum magnitude feedback tap / the 

contribution to the error e[k] is 2 g B . If g fl > 1/2, 

rn<3.x max 

and if I[k] is not ±7, the decision error contribution to 
the symbol estimate pushes the symbol estimate across the 
decision threshold, and an additional error is likely. 

15 When I[k] is ±7, the decision error must have the right 
polarity to push the symbol estimate across the decision 
boundary. When each additional error reaches the same 
feedback tap, a new error is again likely, and a 
potentially long sequence of errors may occur. 

20 The following analysis examines conditions 

under which no error propagation occurs. The following 
simplifying assumptions can be made: the equalizer is 
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unbiased ( gph^^ =1); there is negligible residual ISI 
g T F H RF I F [k] and g£H 

rp2^ P2^k] are negligible) ; and, there is 

negligible noise contribution (g^y R [k] is negligible). 

In other words, the focus of the analysis is 
placed only on the decision errors in equation (20). 
With the above assumptions, equation (20) can be re- 
written according to the following equation: 

f~ \ Nfb 

e[k] = g T B (l[k - 1] - I P1 [k]) = £ g B [«W - 1] - IE* - 1]) (21) 

m = l 

When using a standard slicer as the decision device 22, 
the probability of error is zero when \e[k]\ < 1. This 
condition can be assured if the following inequality is 
true : 

N FB 

2 \g B l*pVc - 1] - I[/r - 1]| < 1 (22) 

m = l 

If it is also assumed that the probability of 
making a decision error, which is not a nearest neighbor 
of the correct symbol, is negligible, then all decision 
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errors have magnitude 2. That is, \T[k] - J[Jc]| < 2 for all 
k. Therefore, if 

IISbIL = ZI^bH < - (23) 

/n = l Z 

5 

then the probability of a decision error propagating a 
new error (under these assumptions) is zero. ||Sb|Ii ^ s the 
1-norm. 

The analysis above indicates how error 
10 propagation may be eliminated if the feedback tap vector 
g B satisfies a constraint on its 1-norm [(g^)^ . The 

discussion below describes how to calculate new equalizer 
coefficients by directly constraining the 1-norm of the 
feedback tap vector. The equalizer coefficients are 
15 calculated by solving the following minimization problem: 

minimize MSE (g F , g B ) (24a) 
subject to ||g F || 1 < M (24b) 

20 where M is the constraint value and the norm is the 1- 

norm. The minimization is performed over both g F and g B . 
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When calculating tap weights that minimize the 
mean squared error MSE subject to an inequality 
constraint such as equation (24b) , the following 
Lagrangian function can be used: 



Ug F , g B , A ) = MSE(g F , g B ) + 2X (flg^ - m), X > 0 



= MSE(g F , g B ) + 2X 



£|g B H - m 



X > 0 



Then, this expression is differentiated with respect to 
10 g F , gB/ and X , the three derivatives are each set equal 
to zero, and the three equations are solved 
simultaneously for g F , gB/ and X . Iterative methods for 
solving this system of equations are well known. 
However, for this calculation of tap weights by 
15 minimizing the mean squared error MSE subject to the 
inequality constraint of equation (24b) to work, the 
inequality constraint function must be dif f erentiable 
(i.e., its gradient must exist). 

However, the drawback in using the 1-norm as a 
20 constraint is that its gradient does not exist when any 
feedback tap is zero. To avoid this situation, the 
constraint is modified as indicated below so that an 
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analytical gradient exists at all points. The modified 
constraint is given by the following inequality: 



M FB 
m = 0 



where the function f (x) is given by the following 
equation : 



f{x) = \ ^ ^ 2 ^4 (26) 



The three constants fo, f2, and f 4 are chosen so that f (x) 
is continuously dif f erentiable . Under these conditions, 
f 0 is arbitrarily chosen to yield a desired value of f (x) 
at x = 0 between 0 and y/2, inclusive. The constants f 2 
15 and f 4 are given by the following equations: 

3 _ f, 



f >=Ty- 2 ^ (27) 

f4 = vi f °- y i) (28) 



-21- 



Attorney 



Docket 
7229 



This modified constraint function is an approximation to 
the 1-norm constraint function and is equal to the 1-norm 
constraint function over most of the range of possible 
values of g B [m]. The region where this modified 
5 constraint function differs from the 1-norm constraint 
function is in the neighborhood of points close to where 
the 1-norm constraint function is not dif f erentiable . 

The difference between this modified constraint 
function and the 1-norm constraint function depends on 
10 the choice of f 0 and y . This difference is shown in 

Figure 2 where f (x) / y is compared to |x| / X as fo/y is 

stepped from zero to 1/2 in increments of 0.05. The 
dashed line in Figure 2 represents the unmodified 
constraint function and the solid lines in Figure 2 
15 represent modified constraint function for various values 
of the ratio f Q /y. 

The Lagrangian function is modified in 
accordance with this constraint to produce the following 
equation: 

20 

U% F , g B / * > = MSE{g F , g B ) + 2X \g B \ x - M), X > 0 



= MSE(g F , g B ) + 2X 



^ f(g B [m] ) - M 
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where f(gB[m]) is specified in accordance with equations 
(26) -(28) and the constraint function M is set as 
described above. Typical useful values of M are 1, 2, or 
5 3. A typical useful values of y is M/1000. A typical 
value for f 0 is 3y /8. Typical values for f 2 and f 4 are 
then given by equations (27) and (28). 

A decision feedback equalizer 40 that 
implements the present invention is shown in Figure 3. 

10 The decision feedback equalizer 40 includes a feed 
forward filter 42 that includes a plurality of feed 
forward taps, a feedback filter 44 that includes a 
plurality of feedback taps, a decision device 46 such as 
a slicer or trellis decoder, a summer 48, and a tap 

15 weight calculator 50 that sets the tap weights for the 
feed forward filter 42 and the feedback filter 44. 

The tap weight calculator includes a channel 
estimator 52 that estimates the channel impulse of the 
channel through which a signal y is received by the 

20 decision feedback equalizer 40. The feed forward filter 
applies the feed forward taps g F to the received signal y. 
The outputs of the feed forward filter 42 and the 
feedback filter 44 are summed by the summer 48. The 
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output of the summer 4 8 forms the output of the decision 
feedback equalizer 40. The decision device 46 makes 
symbol decisions x for the output of the summer 48 and 
feeds these symbol decisions x to the feedback filter 44 
5 which applies the feedback taps g B to these symbol 
decisions x. 

The channel estimator 52 estimates the channel 
impulse response CIR and the signal-to-noise ratio SNR in 
known fashion based on the received signal y and the 

10 symbol decisions x from the decision device 46. A tap 
weight determiner 54 determines the feed forward tap 
weights g F and the feedback tap weights g B from the 
channel estimate and the SNR estimate based on the 
equations and procedures discussed above. The tap weight 

15 determiner 54 supplies the feed forward tap weights g F to 
the feedforward filter 42 and the feedback tap weights g B 
to the feedback filter 44 for application to the received 
signal y and the symbol decisions x, respectively. 

The unconstrained decision feedback equalizer 

20 in theory provides the lowest mean squared error at the 
output of the decision feedback equalizer, as long as 
there is no error propagation. The constrained decision 
feedback equalizer results in a larger mean squared error 
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at the output of the decision feedback equalizer, but 
reduces the chance of error propagation. 

For the constrained decision feedback 
equalizer, as the constraint value M increases, the mean 
5 squared error is lowered, but the effectiveness of the 
constraint is reduced. As the constraint value M 
approaches infinity, the decision feedback equalizer 
effectively becomes unconstrained. As the constraint 
value M is reduced, the chance of error propagation is 

10 reduced but the mean squared error at the output of the 
decision feedback equalizer is increased. If the 
constraint value M is too small, then the mean squared 
error will become too large. In that case, the 
constraint will actually cause error propagation. So, 

15 the constraint value M needs to be confined within a 
useful range. 

Moreover, the constraint value M may be 
adjusted adaptively by an arrangement shown in Figure 4. 
As shown in Figure 4, a block 60 provides an estimate of 

20 the mean squared error MSE (hereinafter MSE' ) at the 
output of the decision feedback equalizer 46 from L 
samples as given by the following equation: 
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MSE 



i L 




where a is the output of the decision feedback equalizer 
40 (the output at the summer 48) and b is the output of 
5 the decision device 46. A block 62 then determines the 
constraint value M according to the following equation: 



10 where a is a desired proportionality constant. Thus, as 
shown in Figure 5, the constraint value M is 
proportionally increased as the estimated mean squared 
error MSE' increases because increasing the constraint 
value M will tend to reduce the mean squared error MSE. 

15 This upward adjustment of the constraint value M will 

prevent the constraint value M from becoming too small. 
On the other hand, as also shown in Figure 5, when the 
estimated mean squared error MSE' is very low, the 
constraint value M is proportionally decreased. This 

20 downward adjustment of the constraint value M will 

increase the mean squared error MSE somewhat, but will 
further reduce the chance of error propagation. 



M = a MSE' 
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The block 62 can also adjust the constraint 
value M according to the function shown in Figure 6. 
Accordingly, if the estimated mean squared error MSE' is 
less than or equal to a low threshold Th L/ the constraint 
5 value M is given as M = a MSE' . However, if the 

estimated mean squared error MSE' is greater than or 
equal to a high threshold Th H , the constraint value M is 
given as M = J3 MSE' . If the estimated mean squared error 

MSE' is between Th L and Th H , the constraint value M is not 
10 changed. Thus, there is a deadband between Th L and Th H . 
The proportionality constants a and J3 may have any 
desired value and may be equal or not as desired. 

Certain modifications of the present invention 
have been discussed above. Other modifications will 
15 occur to those practicing in the art of the present 

invention. For example, the present invention has been 
particularly described above for 8 VSB digital television 
systems. However, the present invention is applicable to 
other types of systems and to other modulation techniques 
20 such as quadrature amplitude modulation (QAM) . 

Moreover, Figure 5 illustrates that the 
constraint value M is a linear function of the mean 
squared error MSE, and Figure 6 illustrates that the 
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constraint value M is a first linear function of the mean 
squared error MSE below the threshold Th L and is a second 
linear function of the mean squared error MSE above the 
threshold Th H . However, functions other than linear 
functions can be used to adapt the constraint value M. 

Accordingly, the description of the present 
invention is to be construed as illustrative only and is 
for the purpose of teaching those skilled in the art the 
best mode of carrying out the invention. The details may 
be varied substantially without departing from the spirit 
of the invention, and the exclusive use of all 
modifications which are within the scope of the appended 
claims is reserved. 
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